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We apply the three-flavor (Polyakov-)Nambu-Jona-Lasinio model to generate baryons as quark- 
diquark bound states using many-body techniques at finite temperature. All the baryonic states 
belonging to the octet and decuplet flavor representations are generated in the isospin-symmetric 
case. For each state we extract the melting temperature at which the baryon may decay into a 
quark-diquark pair. We seek for an evidence of the strangeness dependence of the baryon melting 
temperature as suggested by the statistical thermal models and supported by lattice-QCD results. 
A clear and robust signal for this claim is found, pointing to a flavor dependence of the hadronic 
deconfinement temperature. 


I. INTRODUCTION 

Experiments at the relativistic heavy-ion collider and 
the large hadron collider (LHC) have shown that a quark- 
gluon plasma (QGP) is produced during the first stages 
of a relativistic heavy-ion collision. The QGP is the 
phase of quantum chromodynamics (QCD) at high tem¬ 
perature/density where quarks and gluons are not color- 
confined into hadrons. From this QGP phase, the pro¬ 
duced fireball undergoes a transition to the hadronic 
phase at a given hadronization temperature. 

At nearly vanishing baryochemical potential the phase 
transition to the hadronic state is known to be a 
crossover (I . Experimentally it is known that at the 
so-called chemical freeze-out temperature, which at van¬ 
ishing chemical potential is close to the hadronization 
temperature, the hadrons are in statistical equilibrium. 
This is the result of a fit of the hadron abundances in the 
framework of a statistical model [3 [Si. This fit deter¬ 
mines the chemical freeze-out temperature and describes 
the multiplicity of almost all nonresonant hadrons with 
an astonishing precision. After chemical freeze out the 
hadrons still interact but the chemical composition of the 
hadron gas remains (almost) unchanged. Results from 
high energetic central Pb+Pb collisions at LHC show 
that the freeze-out temperatures extracted by thermal 
fits g] are close to the crossover temperature predicted 
by lattice-QCD studies 0. 

A natural question to ask is whether the freeze-out con¬ 
ditions depend on the hadron species, i.e., if the chem¬ 
ical freeze-out temperature depends on flavor. Thermal 
fits presented in Ref. g] show a tension when fitting the 
different baryonic species with a common freeze-out tem¬ 
perature, suggesting that the chemical freeze-out temper¬ 
ature for nonstrange baryons is smaller (around 16 MeV) 
than that for strange baryons (21 0]. In Refs. [SHE! it is 
shown that thermal fits with two independent freeze-out 
temperatures (for non-strange and strange hadrons) pro¬ 
vide a better description of the hadronic yields and con¬ 
siderably reduce the y 2 (per degree of freedom) of the fit. 
This holds for a wide range of collision energies between 
y/SNN = 6.27 GeV and y' Snn = 2.76 TeV. In particu¬ 


lar, for Pb+Pb collisions at LHC with i/sjvw = 2.76 TeV 
the difference between the two freeze-out temperatures is 
around 15 MeV 16]. 

If this is the case, one may suggest that the hadroniza¬ 
tion temperature also depends on the strangeness content 
of the hadron. This idea was brought up quite recently by 
the authors of Ref. [5]. In this reference, the strangeness 
dependence of the crossover temperature has been stud¬ 
ied with continuum-extrapolated results of lattice-QCD 
calculations. The conclusion was that the crossover tem¬ 
perature (measured by the maximum of a susceptibil¬ 
ity ratio) is about 15 MeV larger for strange hadrons 
than for those composed by light quarks. This differ¬ 
ence is in surprisingly good accordance with the results 
from statistical-thermal fits of ALICE abundances, even 
if the two physical processes (hadronization and chemical 
freeze-out) are conceptually distinct. 

In this paper we study the flavor dependence of the 
hadronization temperature by using one of the simplest 
effective models for strong interactions. The Nambu- 
Jona-Lasinio (NJL) model is an effective model for low- 
energy QCD where the gluonic fields are integrated out 
and the basic interaction consists of a 4-quark contact 
vertex. Although the gluon dynamics is absent in this 
model, some of the gluonic features can be reproduced 
by the so-called Polyakov-Nambu-Jona-Lasinio (PNJL) 
model. 

This effective model lacks true confinement. However, 
hadrons can be thought as dynamically generated states 
from multiquark rescattering, thus providing a nonper- 
turbative mechanism for an effective confinement. The 
properties of these hadrons (masses and widths) can 
be obtained by solving the Bethe-Salpeter (BS) equa¬ 
tion (for mesons) and the Fadeev equation (for baryons) 
with some approximations. Many approaches have been 
applied in which meson and baryon properties at zero 
temperature have been computed within the NJL/PNJL 
models HDHU. These models can be extended to finite 
temperatures and densities. Such an extension allows for 
calculating the “Mott temperature”, the temperature at 
which hadrons are not bound anymore, because they can 
melt into a quark and a diquark. 
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Our aim is to find the Mott temperature for sev¬ 
eral hadrons within the three-flavor NJL/PNJL mod¬ 
els, and extract conclusions about its dependence on the 
strangeness content of the hadrons. In Sec. [II] we intro¬ 
duce the NJL and PNJL Lagrangians and provide a short 
remainder on how a meson can be effectively described 
as a bound state of a quark-antiquark pair. In Sec. m we 
use the Bethe-Salpeter equation for two quarks to gener¬ 
ate diquarks and extract their properties as a function of 
the temperature. In Sec. |IV| we apply the Fadeev equa¬ 
tion to generate baryons as bound states of quarks and 
diquarks. We will present the baryon masses as a func¬ 
tion of temperature up to the melting (or Mott) temper¬ 
ature. We consider all baryons belonging to the octet 
and decuplet flavor representations in the isospin limit 
(up and down quarks with equal masses). Finally, in 
Sec. [V] we present our conclusions, summary and out¬ 
look. We also include six appendices with technical de¬ 
tails, so keep this paper as much self-contained as pos¬ 
sible: In App. [A] we detail the Fierz transformation for 
the NJL Lagrangian and discuss different versions found 
in the literature. Next, in App. [B]we shortly define the 
thermal functions needed in the calculation of the quark 
condensate and the meson/diquark polarization function. 
Appendix[C]is devoted to the simplification of the Bethe- 
Salpeter equation for diquarks and provide the flavor ma¬ 
trices for the different diquark sectors. In App. [D] we 
reduce the Fadeev equation for baryons to a tractable 
form using the “static approximation”. In App. [E] we 
present the projectors on the physical baryon states be¬ 
longing to the octet and decuplet flavor representations. 
Finally, App. [F] is devoted to the simplification of the 
quark-diquark polarization function, providing useful ex¬ 
pressions. 


II. NAMBU—JONA-LASINIO MODEL FOR 
QUARKS 


quarks. In this way, the gauge fields disappear from the 
theory and all what remains is an effective coupling g 
between quarks, related to the strong coupling constant 
and to the gluon mass. 

For three flavors, the resulting effective interaction La¬ 
grangian reads 

c = E(^ , (1) 

a' ij 

where i,j = = 3 are flavor indices and a 1 = 

1 — 1 are color indices with T a being the color 

generators, which for N c = 3 they are represented by the 
Gell-Mann matrices, 

T a ' = \ a ' , (2) 

with the standard normalization 

tr c (T a ' T b ') = 26 a b ' , (3) 


where tr c denote the trace in color space. 

After performing a Fierz transformation maun this 
Lagrangian can be reexpressed in a convenient way to 
describe the qq , qq and qq scattering. Using the Fierz 
transformation described in App. [X] we first obtain the 
qq sector of the effective theory, which will allow us to 
describe mesons. For instance, the pseudoscalar sector 
of the interacting Lagrangian (A4) reads (all repeated 
indices are to be summed) 


£q q = G (ipi i 7 5 t“- *75 rfa ipi) 


(4) 


where a = l,...,iVj — 1 and G is a coupling constant, 
proportional to the original g in Eq. 0. In this work, we 
will take G as a free parameter to be fixed by comparing 
our results with the experimental hadron spectrum. The 
flavor generators r“ follow the normalization 

tr (r a r b ) = 2 6 ab . (5) 


A. Effective Lagrangian 

The NJL model TUI J 1 EL HE describes the low- 
energy interactions of quarks by a four-fermion contact 
vertex. Being an effective model of QCD it respects 
the symmetries of the underlying theory, in particular 
the Uy( 1) 58 SUv(Nf) 58 SUa{Nj) global symmetries of 
the massless QCD Lagrangian. The Uy{ 1) symmetry 
leads to the baryon number conservation, while the chi¬ 
ral symmetry SUv(Nf)®SUA(Nf) is spontaneously bro¬ 
ken down to SUv(Nf) at low temperatures. The £/a(1) 
symmetry is broken by the axial anomaly. 

The form of the NJL Lagrangian is motivated by QCD. 
Consider the quark-antiquark current-current scattering 
in the t channel, mediated by a dressed gluon. The gluon 
propagator gives a factor 1 /(t — nig), with m g being the 
dressed gluon mass. In the limit of low momentum trans¬ 
fer one can neglect t in comparison with the gluon mass, 
thus producing an effective contact interaction among 


For Nf = 3 they can be represented by the Gell-Mann 
matrices. 

The axial anomaly is responsible for the Ua( 1) break¬ 
ing and gives rise to the observed g — rj mass splitting. 
To account for this effect in the our model, we comple¬ 
ment the NJL qg-Lagrangian in Eq. Q with the’t Hooft 
Lagrangian: 

C H = -H det [ijji (I - 75 ) 1 / 7 ] - H det [ 1/7 (I + 75 ) 1 / 7 -] , 

*3 *3 

(6) 

where H is an additional unknown coupling and I is 
the identity matrix in Dirac space. For Nf = 3 this 
Lagrangian represents a six-point fermion interaction, 
which is effectively projected onto a four-fermion inter¬ 
action by using the mean-field approximation [ 10 ) Ho¬ 
using the same approximation, the quark masses obey 
the gap equation 

nii = m i0 -4G(‘tp i 'ip i )+2H('ip j 'ijj j )('ip k il) k ) , j,k^ i\j ^ k 

( 7 ) 
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with rriio being the bare quark mass for flavor i, and the 
quark condensate defined as 

(Aipi) = ~iN c tr 7 J Si(k ) , ( 8 ) 

with the trace acting in Dirac space. We represent by Si 
the dressed quark propagator, 


The final expression for the quark condensate in the NJL 
model is shown in App. [B| 


Euclidean space (with A M = g s A£T a ). The order param¬ 
eter of the deconfinement phase transition (in the absence 
of quarks) is chosen to be the Polyakov loop, <!>, which is 
the thermal expectation value 

$ = p tr c (T) , (13) 

where the trace in taken in color space. To account for 
deconfinement effects via the Polyakov loop, an effective 
potential ZJ(<I>, <E>, T), is added to the effective NJL La- 
grangian £ —► £ —U. U is a function of the Polyakov 
loop and its complex conjugate, which are taken to be 
homogeneous fields. The form of the effective potential 
is inspired by the Z 3 center symmetry m 


B. Medium effects 

To calculate the hadron properties at finite tempera¬ 
ture, we use the imaginary time formalism with the pre¬ 
scription 


/ 


d 4 k 

(2^ 



( 10 ) 


W(T,$,<I) 


rp 4 


Mp-QQ _ h ($3 + $ 3 ) + hi y 


(14) 


with 


T 

b 2 {T) = a 0 + a +a 2 



(15) 


with T the temperature and k° —» iui n the fermionic Mat- 
subara frequencies iuj n = ? 7 rT( 2 ro + 1). 

To account for the finite baryonic density we can in¬ 
troduce a quark chemical potential by adding to the La- 
grangian the term 

£n = ' 52 '*Pi Vij'yoi’j , ( 11 ) 

ij 

where Hij = diag (/z„, fJ.d, l^s) contains the quark chem¬ 
ical potentials (which can be alternatively expressed in 
terms of the baryon, charge, and strangeness chemical 
potentials, HQ, gs)- In this work we will restrict our¬ 
selves to a vanishing chemical potential g u = = g s = 

0 . 


The parameters do, tzi, a 2l 03 , 63 , 64 and To are fitted 
from the pure-gauge lattice-QCD equation of state at 
zero chemical potential m ■ The numerical values of our 
parameters are given in Table [TJ Following the reasoning 
of Ref. P2] we have considered the running of To with the 
number of flavors. As a consequence, we have modified 
the original parameter To = 270 MeV for the Yang-Mills 
case (Nf = 0) to a value of T 0 = 190 MeV for our case 
(Nf = 2 + 1 ). 

This model is called the Polyakov-Nambu-Jona- 
Lasinio model and has been widely used in similar studies 
as ours, e.g., for QCD thermodynamics EZ! or generation 
of bound states [HI]. 

The PNJL grand-canonical potential reads 


C. Polyakov—NJL model 

In the NJL Lagrangian, the gluon fields have been in¬ 
tegrated out of the fundamental theory. However, one 
can still introduce a source of gluonic effects through the 
Polyakov line 


T(x) = V exp 



drAi{T, x) 


( 12 ) 


ttpNjLi®, rm, T) = U{T , $, $) + 2 

i 

d 3 k 


i i 

d 3 k 


r Ej 


2r E 


(2tt) 3 

+ tr c log (1 + L^e~ Ei ^ T 


J (2tt) 3 
i 

tr c log (l + Te- E ‘/ T ) 


( 1 + jL t e -^/T) 


(16) 


where /3 = 1/T, V is the path-ordering operator and with = \Jk 2 + to?. Using the mean-field approxima- 

A 4 = iA°, the temporal component of the gluon field in tion one has [T9| 
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tr c log^l + Le Ei / T 'j 
tr c log(l + lSe~ E */ T ) 


log (l + 3(4) + <5>e- Ei/T )e- Ei/T + e~ 3Ei/T ^ , (17) 

log (l + 3(4> + <S>e- Ei / T )e~ Ei / T + e~ 3Ei/T ^J . (18) 


We minimize the grand-canonical potential with respect 
to the order parameters: $, l>, ( 'ipi'ipi }, 

dflpNJL dVtppijL _ dtlpNJL 

5 $ ’ di- ’ d(^> ' 

The last equation provides the expression for the quark 
condensate in the PNJL model. 


addition, we also show the Polyakov loop as a function of 
T. These quantities serve as order parameters of the chi¬ 
ral and deconfinement phase transitions. The transition 
temperature can be defined as the point at which the sus¬ 
ceptibility (derivative of the order parameter) peaks. In 
the right panel of Fig. [I] we plot the dimensionless chiral 
susceptibilities defined as 


{ipiipi} = -2 N c 

with 


d 3 k 
( 2 tt ) 3 K, 


r [1 -ft(Ei) -/*(£?i)] 


( 20 ) 


($ + 2 ^e~ Ei/T )e~ Ei/T + e ~ 3Ei/T 
1 + 3(4> + ^e~ Ei / T )e~ Ei / T + e _3Ei / T 
(4> + 24>e _Bi ' /T )e _jEi / T + e~ 3Ei / T 
1 + 3(4> + 4>e _ - Ei / T )e _ - Bi / T + e _3Ei / T 

For Nf = 3, the gap relation |7]) and the first two 
equations in (19) form a system of five coupled equations. 


ftiEi) 

U(Ei) 


The system needs to be solved numerically to obtain the 
value of the Polyakov loop (and its conjugate) and the 
quark masses. In the isospin limit, two equations are 
degenerate giving m u = m d (and (V’u'f/’u) = (4 >d^d ))• In 
addition, at vanishing chemical potential one has 4> = $, 
which is evident from our equations. 

In the left panel of Fig. [T] we show the temperature 
dependence of the light and strange quark condensates, 
(VkdAu) and (4>s' l Ps} for both NJL and PNJL models. In 


_ 1 d(ip u ip u ) _ _ 1 d^^s) 

^( 4 ’ulp-u.) — J12 dT > X{$ s ljl s ) — rp2 

The maximum of the susceptibility will indicate the 
chiral transition temperature. They read 7)^= 
246/262 MeV (NJL/PNJL) and T { ^ s) = 238/255 MeV 
(NJL/PNJL). We also show the deconfinement suscepti¬ 
bility defined as 


_ ^ 
» = T if 


(24) 


whose maximum gives the approximate position of the 
deconfinement transition temperature = 181 MeV. 

D. Mesons as bound states of qq 


To obtain the meson propagator one must solve the BS 
equation for the quark-antiquark scattering amplitude i+ 
j —> m + n (latin subindices will denote quark flavor 
and barred indices antiquark flavor) in the random-phase 
approximation (RPA): 


rpab / 2 \ _ K~ a P I A f ^ ^ fC ac 


where a, b denotes the meson flavor channel. The kernel 
K, reads 

K?imn = ^2K ab n b nm , ( 26 ) 

with the vertex factors containing color, flavor and spin 
matrices 

= (icolor ® rg <g> r) , (27) 

as well as a combinatorial factor of 2. The Dirac structure 
whose indices we have omitted in our BS equation can 


Sp (k + f) Sg (k - |) T*l mn (p 2 ) , (25) 


be chosen to be F = { 1 , * 75 , 7 ^, 757 ^} for scalar, pseu¬ 
doscalar, vector, and axial-vector mesons, respectively. 


The addition of the ’t Hooft term to the NJL La- 
grangian breaks flavor symmetry at the mean-field level 
of the coupling constants. The resulting flavor-dependent 
couplings in the qq sector K ab are combinations of the 
G and H couplings m ■ In the pseudoscalar sector, the 
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FIG. 1: Left panel: light and strange quark condensates (for both NJL and PNJL models) and Polyakov loop as a function of 
the temperature. In the legend, u and s stand for (ipulpy,}(T)/(tpuipu}^ = 0) and (ip s ip a )(T)/{ip s ip s )(T = 0), respectively. Right 
panel: Chiral and deconfinement susceptibilities as a function of temperature. Their maxima show the chiral and deconfinement 
transition temperatures. In this panel, the labels u,d, and $ stand, respectively, for X(ip u ip u )(T),X(,p ll ^ s ){T) and y$(T). 


nonzero couplings read mm 

K 00 = G + y ((iputpu) + (V’dV’d) + , (28) 

K 11 = K 22 = K 33 = G - y , (29) 

A 44 = A 55 = G- , (30) 

K 66 = I< 77 = G — y (4> u 1pu) , (31) 

K S8 = G - y + 2(V^d> - , (32) 

K 03 = K 30 = JjL_ - (tp d ip d )) , (33) 

K os = JsT 80 =-^ (( 1 /i„i/’«) + (t/idV’d) - 2(VS s ^ s ))(34) 

A ' 38 = A ' 83 = -^= {{‘tpui’u) - (Edit’d)) • (35) 


Notice that they are diagonal in flavor space except for 
the (0-3-8) subsystem. These non-diagonal couplings will 
eventually bring a n 0 — rf — rf mixing m, which should 
be solved in the coupled-channel basis. In the isospin 
limit ( m u = mf) the n° is decoupled from the system, 


but mixing is still present in the if — rj s subspace. In a 
diagonal basis, this fact accounts for the p — rf mixing, 
providing the mass splitting between these two states. 
Note that in the absence of the’t Hooft term, the mixing 
disappears. 

Introducing the function t ab (p 2 ) 

T£mn(P 2 ) = % t ab (P 2 Mm , (36) 

the solution of the BS equation is a matrix in flavor space 


*“V) = 


2 A 


1 - 2A7I(p 2 ) 

where the polarization function II ab (p 2 ) is defined as 


(37) 


n a V) = * 


d A k 

( 2 *)< ^ 




r< ( fc + f) 




(38) 


In App. [B] we provide a simplified expression for this 
function at finite temperature. 

The poles of t ab (p 2 ) represent the mesonic states, which 
are bound states of the qq scattering. One can perform 
a Taylor expansion of the function t~ 1 ’ ab (p 2 ) around the 
pole p 2 = m 2 M , 


t_1 ' a6(p2) " t " 1 ’ a6(W ^ ) + d± ~dp ^ 1 „ 2 , (/ - ™m) + ■ ■ ■ - - 2mM dp 


1 dU ab {p 2 ) 


C P 2 - m 2 M ) , (39) 


r 


where we have used that t 1,ab (m 2 M ) = 0 at the pole we show that t ab (p 2 ) can be identified with the meson 
position. Defining the effective coupling propagator 


SM->OI7 = 


2 niM 


dn° b ( P 2 ) 


dp 

P 2 = m M 


(40) 


t ab {p 2 ) = 


9 M—>qq 
p 2 - m 2 M 


(41) 
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Parameter 

m q o 

m s0 

A 

G 

H 

Gdiq 

Gdiq.v 

Value 

5.5 MeV 

134 MeV 

569 MeV 

2.3/A 2 

11/A 5 

1.56 G 

-0.639 Gdiq 

Parameter 

do 

ai 

a 2 

03 

63 

64 

T 0 

Value 

6.75 

-1.95 

2.625 

-7.44 

0.75 

7.5 

190 MeV 


TABLE I: Parameters of the NJL and PNJL model used in this study. In the isospin limit we have m q o = m u o = rrido- 


Therefore, the equation, 

1 - 2K ab U a \p 2 = m 2 M ) = 0 , 


(42) 


gives the meson mass tom in the appropriate flavor chan¬ 
nel (and spin channel by selecting T). 

If the generated state has a itim larger than the sum 
of quarks masses, it is possible for this meson to decay 
into a quark-antiquark pair. In this case the polarization 


function (38) becomes complex and the pole acquires an 
imaginary part. Considering the variable p 2 as complex 
one can identify the mass and the decay width with the 


Meson 

7T 

K 

V 

v' 

P 

K* 

U) 


Mass at T = 0 

135 

492 

512 

987 

723 

866 

723 

998 

Mass - pdg [2T 

136 

495 

548 

958 

775 

892 

782 

1020 

Tmou - NJL 

267 

271 

245 

0 

234 

250 

234 

372 

Tmou - PNJL 

282 

286 

262 

0 

253 

266 

253 

382 


TABLE II: Masses at T = 0 and Mott temperatures for 
the pseudoscalar and vector mesons in the NJL and PNJL 
models. For the rj meson we find a finite decay width already 
at T = 0. All units are given in MeV. 


Using the parameter set in Table[j]we obtain at T = 0: 
the light-quark condensate (V’uV’u) = —(241.3 MeV) 3 , 
the pion decay constant = 92.2 MeV, the pion mass 
tti-jy = 134.8 MeV, the kaon mass mx = 492.1 MeV, the 
r\ — r( mass splitting of 475.5 MeV, the proton mass 932.0 
MeV and the A baryon mass 1221.4 MeV. 

Our results for pseudoscalar and vector meson masses 
are summarized in Figs. [2] and [3] respectively, where we 
include the results from both the NJL and the PNJL 
models. 

In Table [TT] we present the masses of the pseudoscalar 
and vector mesons at zero temperature, as well as the 
Mott temperatures for all of them, defined as the tem¬ 
perature at which 

(Tm ott) = rUp(T/vf ott.) J - m.q (Tm ott) •> (43) 

with uip and mq the masses of the quark and antiquark 
that compose the meson. 

Notice that for each mesonic state a different Mott 
temperature is found, note also the large difference be¬ 
tween the Mott temperatures within the NJL and the 


real and imaginary parts of the pole position. In this 
way one can obtain the meson masses and decay widths 
as a function of temperature and/or chemical potential. 
A detailed discussion about this procedure is provided at 
the end of App. [Bj 

The parameters we use in this work are partially based 
on the findings of Ref. [23). For the NJL model in the 
isospin limit we have seven parameters. The extension 
to the PNJL model introduces another seven parameters 
(fixed from the thermodynamics of pure-gauge QCD in 
the lattice |X7J). All of them are summarized in Table |!J 


PNJL models. Given the uncertainties of our model one 
could argue that a common Mott temperature may work 
for them all. Alternatively one might claim that more 
precise data are necessary to establish an eventual differ¬ 
ence of the Mott temperature of pions and kaons, for in¬ 
stance. However, note that the cf> meson has a very large 
Mott temperature, not consistent with the rest. This 
fact makes the </ meson impossible to accommodate in a 
global picture of common hadronization conditions. This 
would manifest itself in a larger </>/7r ratio than the ex¬ 
pected for a production at a common temperature, as 
seen experimentally in Ref. [25j . 


III. DIQUARKS 

A second Fierz transformation can convert the original 
NJL interaction into Lagrangian that describes the inter¬ 
action among quarks HU ■ Their bound states, diquarks, 
belong to a nonsinglet color representation and are not 
experimentally observable states. However, they will be 
important for the construction of baryons. 

It is important to fix the different diquark channels we 
consider. In color space, we neglect the sextet represen¬ 
tation (3 C ® 3 C = 3 C © 6 C ) because the members of this 
representation cannot be combined with an additional 
quark to form colorless baryons (in addition, this chan¬ 
nel is known to be repulsive). In flavor space, diquarks 
from both 3 f and 6 f representations can take part in 
the baryon structure, but they should be appropriately 
combined with the spin structure to have a total antisym¬ 
metric wavefunction mm- A summary of the different 
allowed channels is shown in Tabic Mil All the allowed 
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FIG. 2: Pseudoscalar meson masses as a function of the temperature for vanishing chemical potential in the NJL (left panel) 
and PNJL (right panel) models. 




FIG. 3: Vector meson masses as a function of the temperature for vanishing chemical potential in the NJL (left panel) and 
PNJL (right panel) models. 


Color 

Flavor 

J p 

r 

Denomination 

6s 

Not considered here 


3a 

6s 

1+ 

7m 

Axial 

3a 

3a 

0+ 

*75 

Scalar 

3a 

3a 

0" 

I 

Pseudoscalar 

3a 

3a 

1“ 

7m 75 

Vector 


TABLE III: Different diquarks belonging to different sectors. 
T denotes the spin structure associated with the qq vertex. 


combinations can be alternatively obtained by applying 
a Fierz transformation to the original color-current La- 
grangian into the qq sector mi. All the terms emerging 
from the Fierz transformation, exactly match all the dif¬ 
ferent terms shown in Table III] (see further discussion 
and final Lagrangian in App. AF 

As we will see later, only the low-lying diquarks of each 
spin (scalar and axial-vector ones) will be used to form 
baryons. This is so due to the fact that the masses of 


the pseudoscalar and vector diquarks will be higher than 
the experimental baryon masses and already unstable at 
zero temperature. For this reason we only detail here the 
Lagrangian describing the scalar diquark sector 

C qq = G DIQ ^i l5 Cr A T A '^ T )^ T Ci^T A T A '^) , (44) 
and the one for the axial diquark sector 

C qq = G DIQy (^Cr s T A '^ T )(^ T C 7/i r s T A >) , 

(45) 

where Gdiq and Gdiq,v are coupling constants (related 
to the original g but taken here as a free parameters) and 
G = * 7 o 7 2 represents the charge-conjugation operation. 
Finally, we have denoted by t a and t s the antisymmet¬ 
ric and symmetric flavor matrices, respectively; and by 
T a the antisymmetric color matrices. In particular, the 
presence of the latter reflects that the diquarks cannot 
be color singlets. 

The BS equation for the quark-quark scattering in the 
RPA approximation reads 













































(46) 


T ab (rt 2 ) = K ab 4- ' 
ij,mn / ,V/ ij,mn ' ‘ 


d A k 

(2 W 


/C“^S p (fc + |) S‘(] 


- k 


rjicb 

pq,mn 


0p 2 


where S c (p) = C~ 1 S T (—p)C denotes the charge- 
conjugated quark propagator, with T the transposed op¬ 
erator (not to be confused with the temperature). 

Details concern ing the simplification of this equation 
are given in App. O In terms of the function t ab (jp 2 ), 


T%mn(p 2 ) = Wj *°V) Km , (47) 

we can express the solution of the BS equation as 

2G d/q 


t°V) = yz 


2G DIQ IL“»(p 2 ) 


,2^1 ’ 


(48) 


with the quark-quark polarization function of Eq. (C 61 . 


For the spin-1 channels the polarization function con¬ 
tains a transverse and a longitudinal terms 


n“t = n ab r u „ + m b c h 


where we have defined the projectors, 


Tpv 9 flu 


PpPv 

p2 


r — 


PpPu 

p2 


(49) 


(50) 


The solution of the spin-1 diquark masses involves 
the transverse component of the polarization function 
II ab (p 2 ) = \T^n ab v . For axial diquarks this fact di¬ 
rectly comes from the vector current conservation, but 
for vector diquarks this result still holds m- 
The function t ab (p 2 ) reads 



f a V) 


2 Gdiq y 


ab 


1 - 2GdIQ,V n ± (p2) 


(51) 


where the coupling constant Gdiq,v for vector and axial 
diquarks is, in principle, related to Gdiq by the Fierz 
transformation. However, we will take it here as a free 
parameter to be fixed by a fit to the baryon masses. 

Again, the poles of the t ab (p 2 ) functions are identi¬ 
fied with dynamically generated diquarks in their re¬ 
spective flavor and spin channel. The diquark mass, 
is obtained by the solution of the equation 
1 — 2GDiQ^ ab {rn 2 DI Q) = 0 for spin zero diquarks and 

1 — 2Gdiq,v n±( m l)/Q) = 0 for spin one diquarks. 

In Fig. |4]we present our results for scalar and axial di¬ 
quark masses, which will be used to model baryons in the 
next section. Scalar diquarks are represented by square 
brackets [< 71 ( 72 ] and axial diquarks by parenthesis (< 71 ( 72 ). 
From Fig. [4] it is possible to read off the Mott tempera¬ 
ture for the different states (defined as the temperature 
at which the mass of the bound state equals the sum of 
the quark masses). Beyond this temperature, a diquark 
thermal width is generated, which represents the proba¬ 
bility of the diquark to decay into a pair of quarks. 


FIG. 4: Scalar diquark (upper panel) and axial diquark 
(lower panel) masses as a function of the temperature for van¬ 
ishing chemical potential in the NJL model. For comparison, 
the thermal quark masses are also shown, q stands for light 
quark (q = u = d). 


For completeness, we also present our results for the 
pseudoscalar {< 71 ( 72 } and vector < < 71(79 > diquarks in 
Fig. m As we have anticipated, at T = 0 these states 
have a finite thermal decay width which excludes a role 
in forming stable baryons at low temperatures. For this 
reason, we will neglect these states hereafter. 

We now turn to the PNJL model and show the results 
for the scalar and axial diquark masses with the parame¬ 
ter set presented in Table[|] In Fig.[6]we show the diquark 
masses as a function of temperature at zero chemical po¬ 
tential. In this case, the quark masses are more stable 
as a function of the temperature, generating a system¬ 
atically larger Mott temperature in comparison with the 
NJL model. At T = 0 the masses obtained from the 
NJL and PNJL models coincide, providing a consistency 
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FIG. 5: Pseudoscalar diquark (upper panel) and vector di¬ 
quark (lower panel) masses as a function of temperature for 
vanishing chemical potential in the NJL model. Notice that 
the diquark masses are always higher than the combined mass 
of their constituents. Therefore, they are unstable against de¬ 
cay to two quarks. 


Assuming the pole dominance of the diquark propaga¬ 
tion, we can expand the t ab (p 2 ) function around its pole 
to obtain dJ: 


t“V) = - f [qq] T ■ (52) 

P - m DIQ 

Taking the diquark to be at rest (p = 0) the effective 
coupling gf qiq2] ^ qiq2 is defined as 


P[QlQ2]^QlQ2 


2niDiQ 


dn^Po) 

dpo 


Po=m D iQ 


(53) 


For the axial diquarks we find in the pole approximation 


t ab ’^(po, 0 ) 


P(qiq2)^qiq2 
Po ~ m DIQ 



pPpV \ 

m D/Q / 


with 


5 , (9i®)->gi92 an“ b 


—2m 

(po) 


DIQ 


dpo 


P 0 =m O iQ 


(54) 


(55) 


We now present our results for the effective couplings 
of the scalar and axial diquarks as a function of the tem¬ 
perature at vanishing chemical potential. The NJL re¬ 
sults are given in Fig. [7] and those for the PNJL model 
in Fig. [8} Notice that the Mott temperature clearly coin¬ 
cides in these plots with the value at which the effective 
coupling vanishes. 


IV. QUARK-DIQUARK BOUND STATES 


Diquark 

NJL Tmou (MeV) 

PNJL Tmou (MeV) 

M 

256 

272 

[9 s ] 

273 

287 

( 99 ) 

212 

234 

(9 s ) 

233 

251 

(ss) 

307 

319 


TABLE IV: Mott (or melting) temperature for the different 
scalar [giffc] and vector axial ( 5152 ) diquarks in the NJL as 
well as in the PNJL model. 


check. 

In Table |IV| we present a summary of Mott tempera¬ 
tures (with precision of 1 MeV) for the different states in 
the NJL and in the PNJL model. The PNJL model shows 
higher melting temperatures in all cases. From this table 
one already sees a hint for the flavor dependence of the 
hadronization (Mott) temperature. This temperature in¬ 
creases with the strangeness content of the cliquark. 


In the last two sections we have explored the con¬ 
cept of “hadronization” as a dynamical generation of 
bound states from quark and antiquark scattering. This 
idea -quite successful for the generation of mesons and 
diquarks- can be pushed forward to describe baryons as 
bound states of a quark and a diquark. For Nf = 3 we 
will consider both the octet and decuplet flavor repre¬ 
sentations of baryons. Scalar diquarks (those belonging 
to flavor 3 representation) and axial diquarks ( 6 ) will be 
used to build up the baryon octet and decuplet states, 
respectively, according to the decomposition, 

3 © (3 © 6 ) = (1 © 8 ) © (8 © 10 ) . ( 56 ) 

Notice that the members of the octet also contain 
nonzero contribution from the axial diquarks. However, 
previous results at zero temperature have shown that the 
contribution is almost negligible [T2] , For simplicity, we 
will neglect the axial diquark contribution to the baryon 
octet. 

The starting point to describe baryons is the Fadeev 
equation mm ■ 
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FIG. 6: Scalar diquark (left panel) and axial diquark (right panel) masses as a function of temperature for vanishing chemical 
potential in the PNJL model. 



temperature (MeV) temperature (MeV) 


FIG. 7: Effective diquark-quark-quark coupling, defined in Eq. (531 for the scalar diquarks (left panel) and in Eq. (551 for the 
axial diquark (right panel) as a function of temperature in the NJL model. 



temperature (MeV) temperature (MeV) 


FIG. 8: Same as Fig. [7] but for the PNJL model. 
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Xj’ a (P 2 ,q ) 


dk jk.a.6 


where the baryon wave function is denoted by X Ba and 
it carries a quark index (j). diquark index (j), and a 
possible spin index a. 

The kernel reads m 


l)^\p\q,k) = Qll{P\q)Zf a \ q ,-k) , (58) 

with a first term which accounts for the free quark and 
diquark propagators (see right panel of Fig. [9]) 

Qll{P 2 ,q) = S k (P/ 2 + q) itf{P/ 2 - q) (59) 
and a second term 

Z%' a \q, k) = fijf Sti-q - k ) , (60) 


which represents an interaction with an exchanged quark 
(displayed in the left panel of Fig. [9]). 

We do not attempt here to justify the form of the 
Fadeev equation ( |57| ) and we refer the reader to the orig¬ 
inal papers m 126] to learn the rigorous derivation and 
know more details. 

N ever theless, we can provide a simple motivation for 
Eq. {571: If we denote by G(P 2 ) the full baryon propa¬ 
gator, one can form a Dyson equation by taking Q as the 
leading order approximation (free propagation), and then 
considering Z as the elementary interaction (see Fig. 10). 
The Dyson equation can be symbolically written as 


G = g + gzG , 


(61) 


whose solution reads 


G = 


i -gz 


(62) 


(P' 2 ,q,k)X^(P\k) 


= 0 


P 2 =M | 


(57) 


The baryon masses are now extracted as the poles of the 
baryon propagator, so one needs to solve G ~ 1 X(P 2 = 
Mg) = 0, where A' is the baryon wavefunction. Explic¬ 
itly, 


(1 - gZ)X(P 2 = M%) = 0 , (63) 


at P 2 = Mg, which is a simplified version of the more 
complete Eq. (57). 


The technical simplification of the Fadeev equation is 
given in App. [Dj There, we apply the “static approxi¬ 
mation” for the exchanged quark m- This amounts to 
neglecting the exchanged momentum in Z with respect 
to the quark mass mi. This approximation has been in¬ 
tensively used in other works resulting a very good ap¬ 
proximation (the estimated uncertainty is around 5% as 
quoted in [T5]). 

As shown in App. [D] the Fadeev equation can be re- 
casted into a Dirac-like equation for the baryon wave- 
function evaluated at the baryon mass. For a particular 
baryon channel BB' (we denote by B : B' the physical 
baryon in the initial and final state), the equation to be 
solved reads 


g af3 S BB' _ M BB',a/3(p2) 


P 2 =m 2 b 


= o, 


(64) 


with the matrix (in both flavor and spin spaces) M BB 
introduced in Eq. (|D7[), 


M BB',oq 3(p) = A 
mi 


pt ,b v b> 

( 2 ?r ) 4 H kk 


4 4 r'A 




(65) 


For the members of the baryon octet, M BB simplifies to 
(|D9) 


M 


BB' 


( p ) = — V 


t ,Bq-}B' 


mi 


T>P T k n-j TT 
' kk T jl T lk 11 


kk 


(P) 


( 66 ) 


where the flavor matrices r are given on the top of Ta¬ 
ble |VI[ the projection matrices V are given in App. [EJ 


and the quark-diquark polarization function is defined as 

n feS (P) = - J SxgS k (P ~ 9) %(<?) ■ (67) 

For the baryon decuplet the matrix M BB ,a/3 reads 

pTl) 


M BB'^ = _ pbBpB' ^ n a| ( p) , (68 ) 
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FIG. 9: Left panel: Effective coupling in the Fadeev equation which contains a quark exchange. Right panel: Two-particle 
(quark+diquark) propagator in the Fadeev kernel. 



FIG. 10: Dyson equation for the baryon propagator G = Q + QZ G. 


where the flavor matrices for axial diquarks are given 
in the bottom of Table |VI| and the projection matrices 
are also given in App. 1 The quark-diquark polarization 
function is defined as 


IT 




d 4 q 

( 2 ^ 


Sk(P - q) {q) , (69) 


where its final expression is given in App. [F] 

In su mm ary, the baryon masses are obtained by solv¬ 
ing Eq. (64) with the matrix M BB defined in (66) for the 


members of the baryon octet and in (68) for the mem¬ 


bers of the decuplet. The results for the masses and the 
extraction of the melting temperature for the different 
states are given in the next section. 


V. RESULTS AND CONCLUSIONS 


Our results for the baryon masses at finite temperature 
in both the octet and decuplet representations are shown 
in Fig.[ll]for the NJL model. The results using the PNJL 
model are presented in Fig. |T2] We summarized all the 
baryon masses in vacuum (T = 0) in Table [Vj 

From this example, the axial diquark case is straight¬ 
forward. The only differences in the polarization function 
are the numerical prefactor and the explicit appearance 


of Dirac indices in (67). In addition, note that the dif¬ 
ference in the global sign cancels with the different sign 
in the diquark propagators [cf. Eqs. (52 and ([54])). Fi¬ 
nally, we have neglected the term proportional to p^p" 
in Eq. (54). 


The masses are given as a function of temperature for 
vanishing chemical potential. As baryons are considered 
as quark-diquark bound states, the definition of the de¬ 
confinement temperature is slightly modified to take into 
account the possibility of diquark melting. The “decon¬ 
finement” temperature may be a genuine Mott temper¬ 
ature Tmou (baryon), i.e., when the baryon mass equals 
the sum of quark and diquark masses. Nevertheless, it 
can also happen that the diquark melts at a lower tem¬ 
perature Tmou (diquark) < Tm 0 u (baryon). Therefore, 
we define the baryon deconfinement temperature as the 
minimum of the two: 


T d = mm{T Mott (baryon), T Mott (diquark)} , (70) 


excluding the possibility that the three-quark system is 
stable when the diquark becomes unstable. 

Baryon masses are plotted up to this temperature, 
which is summarized in Table [V| for all the baryon states. 
Beyond Td the baryon mass becomes complex and, in 
principle, one should perform an analytical continuation 
of the Matsubara frequencies to nonreal energies. This is 
beyond the scope of this work. 

Comparison with previous results mi mi] show a 
similar trend for all the masses. In particular the sur¬ 
prisingly low melting temperature for E baryon is also 
captured by 29) ](£(T where, in addition, the melting tem¬ 
perature of proton is slightly larger than the A, as op¬ 
posed to our case. In Refs. [21, 31J the decuplet is shown 
for the first time. However, i n th ese works the Dirac 
structure of the Dirac equation (D7) is simplified by tak¬ 


ing the trace in the quark-diquark polarization function 
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FIG. 11: Baryon masses as a function of temperature for 
vanishing chemical potential in the NJL model. 


Baryon 

P 

A 

E 

s 

A 

E* 

s* 

fi 

Mass at T = 0 

932 

1078 

1152 

1269 

1221 

1366 

1512 

1658 

Mass - pdg (27 

939 

1116 

1193 

1318 

1232 

1383 

1533 

1672 

TMott - NJL 

234 

252 

156 

272 

200 

211 

219 

275 

Tmou - PNJL 

254 

269 

195 

287 

223 

231 

239 

288 


TABLE V: Masses at T = 0 and Mott temperature for the 
different baryons in the octet and decuplet representations for 
both NJL and PNJL models. All units are given in MeV. 


(with no further justification for this step). In addition, 
the transverse nature of the axial diquark propagator is 
omitted. In spite of this fact, the trend for the baryon 
octet coincides with our results, also in the larger mass 
of the A with respect to the proton. The decuplet sector 
is also similar to ours. 

In the decuplet sector we find a hierarchy based on 
the strangeness content, i.e. the T c increases with the 
number of constituent strange quarks of the baryon. In 
this respect the first baryon that melts is the A, whereas 
the SI has the largest Mott temperature. 

The main conclusion of this work is the evident fla¬ 
vor dependence of the deconfinement temperature of 



FIG. 12: Baryon masses as a function of temperature for 
vanishing chemical potential in the PNJL model. 


baryons. On the qualitative level, our result is quite ro¬ 
bust due to the fact that the relative temperatures remain 
similar for both NJL and PNJL models. On the quan¬ 
titatively side, the Mott temperatures are strongly de¬ 
pendent on the model (NJL or PNJL) used, and in both 
cases they seem to overestimate the standard values for 
the freeze-out temperatures obtained by the thermal fits. 

Statistical models applied to ALICE data predict a 
lower chemical freeze-out temperature for proton in com¬ 
parison with the that for states with multistrangeness (5 
and fl). Within our approximation, our findings agree 
with this claim being the proton temperature 38 MeV (33 
MeV) less than the 5 temperature in the NJL (PNJL) 
model. This fits well to the experimental findings in [J]. 

In addition, we find that the temperature of the 5 and 
Q baryons are surprisingly similar, in accordance with 
the thermal-statistical model result [1]. In our scheme, 
this result is totally nontrivial, because of the fact that 
the internal baryon structure is entirely different for the 
two states, because we have different kinds of constituent 
diquarks. 

In summary, we have applied the NJL and PNJL mod¬ 
els, together with different many-body techniques - 
Bethe-Salpeter and Fadeev equations— to generate di- 
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quarks and baryons, respectively. The parameters of the 
models are fitted to agree with the low-lying states in the 
respective channel ( Gdiq fitted to reproduce the proton 
mass and Gdiq,v to the A baryon mass), being the mass 
of other baryons predictions of the models. 

Extending the method to finite temperature we were 
able to predict the temperature dependence of baryon 
masses for all the physical states belonging to the fla¬ 
vor octet and decuplet representations. We find a strong 
dependence on the melting (or deconfinement) temper¬ 
ature depending on the flavor content of the baryons. 
In a qualitative way, our findings coincide with the sug¬ 
gested results on the statistical thermal model on ALICE 
data (4], and supported by lattice-QCD results by [9). 
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Appendix A: Fierz transformation 


The Fierz transformation allows us to convert the 
original NJL Lagrangian in Eq. 0 — based on the 
color-current interaction- into a Lagrangian where the 
fermion fields are reordered to account for interactions in 
different color, flavor and spin sectors. 

Even if the Fierz transformation is a well-known 
tool [[52] , some differences in its application can be found 
in the context of the NJL model. The Fierz transforma¬ 
tion itself is based on algebraic identities in flavor, color 
and spin spaces. It is therefore unique once Nf and N c 
are fixed. However, there are at least two methods to 
obtain the transformed Lagrangian. 

The first method requires two different Fierz transfor¬ 
mations m to obtain all the mesonic and diquark contri¬ 
butions. A Fierz transformation transforms the original 
color-current interaction to the “exchange sector” con¬ 
taining the qq interaction in both color singlet and color 
octet representations. In color space the transformation 
reads: 


AC—1 


E Tv'j'Tk'v = 2 


N 2 c ~l 

m 


fii'l'fik'j' 


! ^-1 


N c 


a '=1 


rpa rpa 

1 i'l' 1 k'j' • 

(Al) 


To study physical mesons, one considers the first term, 
while the second term is simply neglected. A second Fierz 
transformation is applied to the color-current interaction 
to generate the qq sector in both color antitriplet and 
sextet representations: 


AC-1 

E ^ 

a'—l 


rpa rpa 
1 i' j' 1 k'l' 


Nc-l T<S' T’S' N c + 1 rpA' rpA’ 
nr 2 -^/ 1 k ' 1 1 'i’ AT 2^ li,k ' ll 'o' 

C 0 / C A/ 


(A2) 

where S' and A! run over the symmetric and antisymmet¬ 
ric members of the color representation, respectively. For 
N c = 3, S' = 1,3,4, 6, 8 and A! = 2,5, 7. The first term 
is neglected as it contains a repulsive diquark interaction 
and is not useful to generate colorless baryons. 


The second method makes use of only one Fierz trans¬ 
formation na on which follows from the principle of 
obtaining only attractive-color interactions. With this 
prescription, the final Lagrangian for both mesons and 
diquarks is obtained by a single Fierz transformation in 
color space: 


AC —1 

E rj-iO,' rriO,' _ c ^ f C ^ Y ^ rriA’ rriA' 

- 2N di ' lld i' k ' ■ 

a' = l c c A' 

_ ( A3) 

In the right-hand side of Eq. (A31, the first term produces 
the Lagrangian for mesons in the color singlet represen¬ 
tation, and the second one the Lagrangian for diquarks 
in the color antitriplet-color representation. No sign of 
any qq interaction in the color octet and of the qq one in 
the color sextet is seen. As claimed, these channels are 
not needed to obtain physical mesons and the baryonic 
states. 


Notice that the numerical factors appearing in front 
of the Fierzed Lagrangian depend on the choice of the 
method. Therefore, they present different values for the 
coupling constants. However, we remind that we are 
treating these couplings as free parameters to be fixed 
by reproducing the hadron masses. For this reason the 
two methods are equivalent, if the coupling constant is 
considered as a free parameter. 

In this work we will use the conventions given in 
Ref. m and generate all possible interactions in the qq 
and in the qq sectors. The terms which we obtain are 
in exact correspondence with the irreducible representa¬ 
tions of the flavor and color group products. 

For the meson sector the Fierz-transformed Lagrangian 
reads [TT] (we suppress the flavor indices i,j): 
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C 


ex 


2 % a Np 9 [ WO 2 

[ i^T a 1p) 2 

Njs- 0 g[$T a 'ii))‘ 

2W~J [ $r“T“>) 


+ (V’ns V') 2 

+ (^i7 5 r a '0) 2 

+ (^?7 5 T a V) 2 

+ (ipi'Y 5 T a T a 'ip) 


K^r^) 2 

\{ip^T a 'ip) 2 
\{^T a T a ' ip) 2 


\{ipYl ^) 2 ] 
|(^75T a V’) 2 ] 
i(^7 5 T“» 2 ] 
l(^7 5 r“T“>) 2 ] , 


(A4) 


where a = 1,..., IV 2 — 1 and a' = 1,..., IV 2 — 1. It is not 
difficult to check that the different terms match with the 
representations spanned by 


duce the physical mesons. In particular, the second term 
is the relevant one for the pseudoscalar mesons (77 K and 

n)- 


(3 ® 3) c ® (3® 3)/ (g) { 1 , * 75 , 7 ^, 757 ^} . (A5) 

The second row (in bold font) will be of interest to pro¬ 


Performing the second Fierz transformation one ob¬ 
tains the diquark sector M: 


C 


QQ 


^2N^9 [ (^i75Cr A T A V r )(V’ T Ci75'r A T A ''i/>) 

\ (^7^75 Cr A T A 'ip T ) (tp T C'y fl - 75 t a T a ' ip) 
^m }9 [ (dpi^CT S T s ' ip T ){ip T Ci^T S T s ' ip) 

| (V’7 M 75C'T S ? lS ''!/i T ) {%p T C')^T S T s ' ip) 


(ipCT A T A ' ip T ){ip T CT A T A ' ip) 
i('^7 At Cr s T A i/i t )(^ t C7 (1 t s T a ip) } 
(PpCT S T s ' ip T )(ip T Ct s T s ' ip) 
\{PP^Ct a T s 'P p T ){ip T C ll ,T A T s 'ip)} , 


(A6) 


where the indices A, A' refer to the antisymmetric mem¬ 
bers of the flavor and color representations, respectively, 
and the indices S , S' to the symmetric elements of the 
flavor and color representations. In this way, the Fierz 
transformation in the qq sector generates the eight pos¬ 
sible terms consistent with the direct product of 


A [ (3 <g> 3) c ® (3 < 8 ) 3)/ <g> {l,*75,7 M ,757 Al } ] , (A7) 


where A denotes the antisymmetrization operator to re¬ 
spect th e Pa uli principle for the exchange of two quarks. 
In Eq. (A 6 ) we highlight the two terms giving rise to 
the scalar and axial diquark discussed in this work when 
forming baryons. 


Appendix B: Quark condensate and meson 
polarization function 


For completeness, we include here some reduced ex¬ 
pressions for the quark condensate in Eq. ([81) and the me¬ 


son polarization function in Eq. (38). For practical rea¬ 


sons we remind here the A and B 0 functions introduced 
in Ref. [521 G2] ■ These functions are convenient for im¬ 
plementing numerically the momentum integrations for 
the thermal averages. 


For vanishing chemical potentials, the function A is 
defined as 


A(rrii, T, A) = 16tt 2 T^ 

n 


d 3 k 1 

(27t) 3 (iu n ) 2 - E 2 


(Bl) 


with the quark energy Ej = k 2 + m 2 and the fermionic 
Matsubara frequencies iui n = i(2n+ 1 )ttT. This function 
naturally appears when computing the quark condensate 
defined in Eq. (| 8 | : 


— m ■ 

(ipiipi) = N c -^-AA(mi,T, A) . (B2) 

Performing the Matsubara summation [55] one obtains 

/ d 3 k 1 

(2^2E“ [1_2nF(i?l)] (B3) 

= 4 f dEi J E 2 - m 2 [n F (Ei) - n F (-Ei)\ , 


where np(Ei) = (e Ei ^ T + 1) 1 is the Fermi-Dirac dis¬ 
tribution function and A^ = \/A 2 + m 2 . For the PNJL 


The function Bq is defined as 


model one simply replaces Eq. |B2| by Eq. 20 
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B 0 (mi,mj,iv m ,p,T,A) = 16n 2 T ^ 

n 


d 3 k 1 1 

(27r) 3 (iw n ) 2 - Ef ( iuj n - iv m ) - .E 2 


(B4) 


with Ei = *Jk 2 + mf and Ej = y/(k — p) 2 + m 2 . It nat¬ 
urally appears in the quark-(anti)quark loop function, 
when computing the diquark (meson) polarization func¬ 
tion. After performing the Matsubara summation, the 
Matsubara frequency iv rn is analytically continued to real 


values of the energy with the prescription iv rn —► p 0 + ie 
For instance, the pion polarization function ( nrii = m u , 
mj = m d ) can be expressed as a combination of the A 
and B 0 functions as 


A ud {po,p,m u ,m d ,T,A) 


{A(m u ,T, A) + A(m d ,T, A) + [{m u - m d ) 2 - pi + p 2 ] B 0 (m u ,m d ,p 0 ,p,T,A)} . 

(B5) 


In the case of a pion at rest, p = 0, one can perform the Matsubara summation and finds 


where 


777,;. »n. 0. T. A) 


(B6) 


5? B^(mj,mj,p 0 ,T, A) 

S BQ(mi,mj,p 0 ,T,A) 


ip. 

Po 
27T 


•V. [ E dEiE 2 - m 2 

J m. 


TIf(±-Ei) 

Ei — Eo 


2 7 r / 

— \J E l - m i n F {±E 0 )<d{A E - E 0 )Q(E 0 - m») 


(B7) 

(B8) 


with So = -(pq + w? ^ rn 2 j )/{2p 0 ). 

Let us briefly discuss some of the prescriptions to define 
the mass and decay width of the mesons and diquarks. 

These generated states are identified with the poles of 
the scattering amplitude t ab (p 2 ), or Eq. ( |42| ). We can dis¬ 
tinguish those states that are generated below and above 
the two-quark mass threshold, that is, those which can¬ 
not decay into a pair of quarks and those in which this 
decay channel is open because their mass is larger than 
the combined mass of the constituents. 

In the first case we talk of “bound states” (the decay 
width is exactly zero). For them, the polarization func¬ 
tion is a real function ( B 0 does not develop any imaginary 
part) and the pole is generated on the real axis of the po~ 
plane (in the first Riemann sheet). The real value of po 
is associated with the mass of the bound state. In the 
second case we denote them as “resonances” (a finite de¬ 
cay width is generated), the polarization function is now 


complex [Eq. (B8 is nonzero) and the pole emerges at a 
complex po (in the second Riemann sheet). The imagi¬ 
nary part of the pole can be related to the decay width 
of the resonance. 


At finite temperature, the variable po = iv m is strictly 
a Matsubara frequency. If a bound state is generated, 
then one can simply make the standard analytical con¬ 
tinuation to real energies iv m —> po +ie and find the value 
of the generated mass by solving 1 — 2I\ ab Il ab (po) = 0, 
with real II ab (po)- However, for a resonant state, one 
must analytically continue the Matsubara frequency to 
complex energies and find the pole in the second Rie¬ 
mann sheet. As this procedure might be cumbersome (in 
particular for cases with several coupled channels, where 
additional Riemann sheets must be considered), approx¬ 
imate methods are used (see, for instance, BSD- 

One introduces the spectral density p(j>o,p) (see [56] , 
for instance) as the imaginary part of the bound 
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state/resonance propagator, 

P ab (Po,p) = --%t ab (p 0 ,p) (B9) 

7 r 

i_ sn ab (p 0 ,p) _ 

7T [(2G)" 1 - U ab (p 0 , p )) 2 + (QTI ab (po, p )] 2 

where the Matsubara frequencies has been analytically 
continued to real energies. Therefore, the spectral func¬ 
tion is a real function of real argument. 

Note that taking p = 0, using the pole approxima¬ 
tion and the optical theorem we can check that the spec¬ 
tral function is proportional to the scattering amplitude 
squared evaluated on the real axis, 

P ab (po, 0) = -- Zt ab (p o ,0) <x |T| 2 (p o ,0) . (BIO) 

7r 

Thus, whenever the scattering amplitude presents a pole, 
this is reflected into the spectral function as a peak. If 
the pole is not far from the real axis (3po •C Jipo), the 
real part of po at the pole position coincides with the 
maximum of the spectral function. Therefore, the mass 
can be defined as the position of the spectral density 
maximum. If, in addition, 3II ab (po) is a smooth function 
of po [36] around the peak, it can be approximated by a 
Lorentzian shape, with a width mm 

r = -sn ab (p 0 ). (Bii) 


However, this approximation might break down if the 
pole position is far away from the real axis, and in partic¬ 
ular, if the many-body equation presents coupled chan¬ 
nels: a broad resonance can be hidden by another pole 
with more strength in this channel, several poles might 
appear very close in the complex plane producing a com¬ 
bined shape of the the spectral function, a new threshold 
opens close to the resonance (Flatte effect) blurring the 
Lorentzian shape, etc. 


As in this paper we are not interested in a precise ex¬ 
traction of the decay widths but only in the temperature 
at which they become nonzero, we use an intermediate 
prescription described in Ref. [22J [34]. In this case one 
makes the analytical continuation to real energies (after 
having performed the Matsubara summation) but con¬ 
siders a complex po i n the factor in front of the Bq func¬ 
tion in Eq. (B5). This prescription provides a complex 
polarization function of complex argument and one has 
direct access to the mass (m = Sftpo) and decay width 
(T = —23p 0 ) of the resonant state. 


Appendix C: Bethe-Salpeter equation for 
quark-quark scattering 


In this approximation, one defines the decay width of the 
resonance as the Lorentzian width T. 


The Bethe-Salpeter equation for the gg-scattering m 
in the scalar channel reads Eq. (461 


r-nO.b / 2 \ _ 1C a ® 

ij.mn\P ) ij.mn ' 


drk 


(27 


— ic ac 

4 *'ij,pq 


S P 0 + f) ^ c (| -k) Tp b mn {p 2 ) 


(Cl) 


We can pull out all the vertex factors by defining the 
diquark propagator t a b 

T£ mn (p 2 ) = ft“.i“V)o b ro , (C2) 

where 

= r a ' <g> t“. <g> re , (C 3 ) 


as there is no flavor singlet in this channel and its effects 
are expected to be much suppressed (there is a small ef¬ 
fect around 4 %, discussed in Ref. |12l). 

In terms of the amplitude t(p 2 ) we can express the 
solution of the BS equation for scalar diquarks, 


f a V) 


2G_d/q 

1 - 2Gp /Q n° b (p2) > 


(C5) 


with T = 7757 ^ for scalar and axial diquarks, respec¬ 
tively; G = * 7 o 72 being the charge conjugation operator. 

The kernel is taken from the Fierzed NJL Lagrangian 
of Eq. ([441 


_ r\a (~)fr 


(C4) 


where the 2 is a combinatorial factor arising when at¬ 
taching the external legs to the vertex and Gdiq must 
be substituted by Gdiq,v f° r axial diquarks. Note that 
we neglect any contribution from the’t Hooft Lagrangian 


with the polarization function 

d 4 k 


n a V) = * 


(2t r) 


4 Tr [^Mk+p^^Sjip^-k)] , 

(C6) 

where the trace is to be taken in color, flavor, and Dirac 
spaces. For axial diquarks, the amplitude involves the 
transverse part of the polarization function [cf. Eq. (49)], 

2Gdiq,v 


t ab (p 2 ) = 


1 - 2G DIQ yUf{p 2 ) 


(C7) 
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For the scalar case one has 


n“V) = tr c{T a 'T b ') i J -0^ tr 7 [* 7 & (k + |) rysCSf (| - fc) C" 1 


(C8) 



TABLE VI: Flavor matrices for the antitriplet and sextet 
flavor representations of the direct product 3/ ® 3/ of 5(7(3). 


The color factor is common to all diquarks in the an¬ 
titriplet color representation, where the generators can 
be taken as m 


(T“ )y k , = ie a 'j' k ' ■ (C9) 

The color factor reads 

tr c (T a T b ) = ~Y^ e a ' k ' l 'e b ' l ' k ' = 25 a ' b ' , (CIO) 


k’l’ 


which means that the color of the diquark does not 
change in the propagation. From now on we will sup¬ 


press the color indices. 

Using the identity, 

CSJ (| - k) C- 1 = (k - f ) , (Cll) 


we can express the polarization function as 


n“V) = 2T?4t J 


d 4 k 

(27T 


4 tr 7 [i”f 5 Si (k) i^Sj (k - p)] 


(C12) 

where we have performed a variable shift k —> k — p/2. 

In flavor space, we choose the representations shown in 
Table VI for the sextet and antitriplet case [12IHH]. For 
both of them, the normalization is tr ( T a r b ) = 28 ab . 

As an example, we calculate the lightest a = [ud] di¬ 
quark. It is easy to see that the polarization function 
is diagonal in flavor, so the only possibility is to have 
b = [ud]. Using r]“ d] = -A?-, 


n M (p 2 ) = 4i J 0A tr 7 [i^Sd (fc) i'ysSu (k - p)] . 

(C13) 

This is the final expression for the polarization function 
that gives rise to the propagator of the [ud] diquark. 

At finite temperature, we introduce fermionic Matsub- 
ara frequencies, 


_ r fp 

n [ “ dl (w m ,p) = -4T ^2 J tr 7 [i^Sd (iu} n ,k)i'y 5 S u (iu n - iu m , k - p)] . 


(C14) 


r 


After performing the Matsubara sum, we analytically 
continue the unsummed Matsubara frequency iv m to the 
real energy p 0 + ie. For practical purposes, this function 


can be reduced in terms of the A, Bq functions defined 
in (22! [34] and detailed in App. |B| 


n M1 (p 0 ,p) 


l 

2tt 2 


{A(m u , T, A) + A(m d , T, A) + [(m- 


m d ) 2 -Po + P 2 ] Bo(m u , m d , p 0 , p, T, A)} . 
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Similar expressions can be found for other diquarks in 
different flavor channels. 


For axial diquarks, the expression for the transverse 
part of the polarization function reads: 


n^OcP) = {A(m u ,T, A) + A(m d , T, A) + [[m u - m d ) 2 - 2 m u m d - pi + p 2 ] B 0 (m u , m d ,p 0l p, T, A)} . 


37T 2 


r 


Appendix D: Reduction of the Fadeev equation 

In this appendix we will give some details of the simpli¬ 
fication of the Fadeev equation in Eq. (57) and its reduc¬ 


tion to a Dirac equation. We follow the same reasoning 
as in Ref. E3, but with a different notation. 


We start by considering Eq. (571 


9 "Oj'k 


S '' 1 "'6 - 


d 4 k 


(2tt) 


- A Lf k ^{P\q,k) 


Xf(P 2 ,q) 


P 2 =Ml 


= 0, 


(Dl) 


r 


where its kernel was defined in Eq. (581. Notice that we 


have denoted the color terms of the equation by primed 
indices, to distinguish them from flavor factors. 

For convenience, we define a new baryon wave function 
Y by integrating over the momentum q na. 


Y k k 0 (P 2 ) = j 


d 4 q 

(2 W 


Xt\P i ,q) • 


(D2) 


To express the Fadeev equation in terms of the new wave- 
function one integrates Eq. (Dl I over q to get 


5 jk S~ fk Y k \P 2 ) - 


d 4 q 


x r“ S k (P/2 + q)itf{P/2-q) X K k p {P\k) 


fc/3. 


d 4 k 




P 2 = n 


= 0 


where we have substituted the kernel given in Eq. (581. 


Notice that it is not possible to express the equation only 
in terms of Y k ^(P 2 ) because there is a quark propaga¬ 
tor that depends on the momentum k. Therefore, this 
equation is nonseparable and can be only solved with 
numerical techniques [28] . 

In the so-called “static approximation” m one ne¬ 
glects the momentum dependence of the quark propaga¬ 
tor by assuming that the dress quark mass is much larger 
than the typical k: 

S ^-q-k) = —rAi -I- (D3) 

- p- mi mi 

This approximation makes the Fadeev equation sep¬ 
arable, and allows for a trivial integration on k. The 
equation is reduced to the simpler form: 


g 


lk >^ k 'S jk 6~ jk 


_j 4 i _ 1 

a y pk' rpj' k j pMp _ 

(2tt) 4 *'* *'*' * ,fc M 


X S k 


it' 


Oi(3 


Y k ?(P 2 ) 


= 0 . 


P 2 =m% 


To obtain the baryon masses of the different states, we 
project this equation onto physical states B and B’. The 
flavor projectors are defined in App. [E] for those states 
belonging to the octet and the decuplet representations. 
Applying these projectors we find 


SjkS 3 * Vffvg = V^fVfi = 5 BB ' . (D4) 

In color space we take the projector onto the singlet state, 
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V",y >r = S^j'/V 3. In the first term of Eq. (D4) one has 


yrr 

ycolor, f 


AT’ T <W<A' n C ^ r,t = o tr - 1 = !. (D5) 


3 3' 


whereas in the second term of Eq. (D4) 



- 2 , 


(D 6 ) 


where we have used Eq. (C9). 


To simplify the notation we can define the matrix 
M bb ' 


M BB '’ a/ 3 (P) = — 

mi 


— — pt.Sp-B' 

(2t r) 4 P fcfe 


-k j p^ip 

r ji r Zfc 1 1 P 




(D7) 


The final equation is expressed as a Dirac-like equation, 

g »fi 6 BB' _ M BB', a ^p2 = M 2j = 0 ) ( D8 ) 

which is a matrix equation in Dirac and flavor spaces. 


For the baryon octet (composed by scalar diquarks) the 
a,j3- Lorentz indices are absent and T^ = ij 5 . Therefore 
one has (after performing the change of variables q —> 

-q + P/2) 


M bb '{P ) = 


_-piT-p/J 

' A A / U 

mi 33 


d 4 q 

(2^ 


Sk(P-q) it- k (q ) 


(D9) 


where the scalar diquark propagator t k is taken from 
Eq. (pi). 


In flavor space this equation is diagonal except for the 
A — E° -1 mixing, where one needs to solve 


/ i _ m s ° s ° —M s ° a -M ^ 1 \ 


det 


-M as ° 1 - m aa -m ai 
y —M IS ° -M ia 1 - M w J 


= 0 . (DIO) 

In the isospin limit ( m u = nid ) one has = 


M s 1 = M as ° = M s A = 0. Thus, the E° decouples 
from the system (DIO), and in this limit its mass is de¬ 
generate with the mass of the E + baryon. However, the 
elements M IA and M AI are nonzero and this produces 
a mixing between the flavor singlet and the A. There¬ 
fore, to obtain the mass of the A we need to solve the 
two-channel equation. 


For the members of the decuplet (composed by axial 
diquarks) one uses T M = 7 ^ to get 


M bb'^{P) EE — V]fV B i rjj 4 


mi 


A 

(2A 


Sk ( P - q) itr P (q) , 


(Dll) 


with the diquark propagator taken from Eq. (54). In this 
case, the Eq. (D 8 ) is diagonal in flavor space. 


Appendix E: Physical baryon projections 

Baryon masses are computed by solving the Fadeev 
equation projected into the different physical states. The 


baryon projectors project the general wave function onto 
the wave functions of specific baryons. In our scheme, 
the baryon wave functions are the direct product of quark 
and diquarks wave functions. In this appendix we provide 
the precise expressions for completeness ESS HQ]. 

In Table |VII| we present the baryon octet and decu¬ 
plet wave functions in terms of the quark-diquark states. 
As in the main text, scalar diquarks are represented by 
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Octet member 

Wave function 

Decuplet member 

Wave function 

V 

| u[ud]) 

A ++ 

| u(uu)) 

n 

| d[ud\) 

A+ 

A=(|d(ww)) + V2\u(ud))) 

A 

^g(|u[ds]) + |d[us]) - 2 |s[ud]» 

A 0 

^(\/ 2 | d(ud)) + | u(dd))) 

E+ 

Hus]) 

A" 

| d(dd)) 

E° 

7^(1 u[ds]) - |d[us])) 

E* 

~^(\s(uu)) + y/2\u(us))) 

E _ 

|d[ds]> 

E*o 

7 s(kM)> + M(ms)) + |w(ds))) 

77 O 

\s[us]) 

E*- 

7 j(|s(dd)) + V2\d(ds))) 

E _ 

\s[ds]) 

^*o 

W(^|s(us)) + | u(ss))) 

Singlet 

Wavefunction 


^j(^ 2 | s(ds)) + | d(ss))) 

I 

^(|u[ds]> + | d[us]) + | s(ud\)) 

Q- 

\s(ss)) 


TABLE VII: Baryon wave functions for all the members of the octet and decuplet flavor representations. 


square brackets and axial diquarks by parentheses. We contribution to the members of baryon octet, 
remind that in this work we neglect the axial-diquark 


These wave functions help us to construct the baryon 


projectors. They satisfy 

VP = V% , (El) 

where B represents the physical baryon state. They are 
orthonormal within the same representation, 

(V^Vfj = 6 BB ' . (E2) 

For the baryon octet they read 

K = ' < E3 > 

n, = . (E4) 

H = H • < E5 > 

P S" = 75 = \/H • <E6) 

*§* = l(A‘:FiA% , (E7) 

J>5 = l(A 6 + a% , (E8) 

- !(**+«%• < e9 > 

Finally, we need the projector for the flavor singlet state, 

n 3 = \ j \ % ■ (eio) 


Note that for the baryon octet, the diquark index i runs 
from 1 to 3 because the scalar diquark belongs to the 


antitriplet representation of 517/(3) x 517/(3). In partic¬ 
ular, the [ds],[us] and [ud] diquarks are represented by 
i = 1,2,3, respectively. 

For the baryon decuplet the projectors read 


V A++ _ 
ij 

<5/1 <5/1 , 

(Ell) 

71 

■3: > 

+ 

II 


(E12) 

II 

O 

[2 1 

Y ^i2^j2 + yf ’ 

(E13) 

idA — 

ij 

<Mj2 , 

(E14) 

_ 

ij 

1 [2 
■^|(5/i<5/3 + Y 2^14^1 , 

(E15) 

\-\*0 

— 

ij 

~^(fii2$j3 + + S-M , 

(E16) 

— 

ij 

+ ^8is6j2 , 

(E17) 

II 

0 

w,:? 

[2 1 

Y 2^14^3 + yf ’ 

(E18) 

v?* 

IJ 

[2 1 
= Y 3^*5^j3 + ’ 

(E19) 

ij 

= <^6^3 ■ 

(E20) 


In this sector, the axial diquarks belong to the sex¬ 
tet representation of 517/(3) x 517/(3). Therefore, the 
index i runs from 1,...,6 representing the diquarks 
(uu), (ud), (dd), (us), (ds) and (ss), respec tivel y. 

All the projectors satisfy explicitly Eqs. (El) and (E21. 
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Appendix F: Quark-Diquark polarization function 

We will detail here the reduction of the quark-scalar 
diquark polarization function in Eq. (67). At finite tem¬ 
perature, the polarization function (671 reads 


n fcfcOl>P) = T Y 1 


d 3 q 


(27r) 


where ivi is a fermionic Matsubara frequency which will 
be analytically continued to real values at the end of the 
calculation. The variable iui n is a bosonic Matsubara fre¬ 
quency appearing in the diquark propagator in the pole 
approximation (521, 


3 S k (wi-ico n , P - q) f £ (zw n ,q) , 

(FI) 


n fefc(*^» p ) 


P 9[qq]^qq 


d 3 q 


(2tt) 3 (iu n + e kq ){iu n - e kq ) (iu„ 


(Pq - ivnho + ni k I _ 

-Po + E kiq )(iu n - Po - E kiq ) 


(F2) 


With e I,g = m DIQ,k + C l 2 aild E k,q = m l + ( P ~ q ) 2 ( n0te 

that k represents the quark flavor, not a momentum). 

The Matsubara summation is performed taking into 
account the four poles using standard techniques [35]. 
We can express the final result in terms of four functions 
(one coming from each pole): 

U k - k (p 0 ,P = 0) = + Jl + J* + Jf) , 

(F3) 


where we have considered the baryon at rest P = 0 and 
performed the analytical continuation ivi —> P 0 + ie. We 
have also defined four J functions, whose real and imag¬ 
inary parts are given by 


II 

tq- 

1 

=c 

2 fb 

= 

7r 


2 Po 

II 

cq- 

£ 

1 

□= 

2Po 

II 

cq- 

6 

7T 

n= 

2P6 


P.V. / dE kA [1 - 2 n F {±E k ^ q )\ (I lm k =F 7o E k , q ) 


J k,q 


E k ,q — -E 4 


(I m k =F loE ± )\J{E ± ) 2 - m 2 k [1 - 2 n F {±E±)\ 0(A B - E±)Q(E± - m k ) 


/'Ar* 


P.V. 


de k [1 + 2 n B (tk)\ [T+ 7 o(e fc T P 0 )] 


~ m m Q ,k 


L DIQ,k 


£k ~ e± 


[Tl m k + 7o(e± T Po)} \]( e ±) 2 - m D IQ p t 1 + 2n s( e ±)] 0 ( A -E* ~ e±)©(e± - m D iQ,k) 


(F4) 

(F5) 

(F6) 

(F7) 


with A e 

= \/A 2 +m£,A B . = A 2 + rn 2 DIQ j_. 

In addi- 

tion, 

„ , m DIQ,k~ m l- P 0 

Ej-i- = it--- , 

2 Po 

(F8) 


m DIQ,k ~ m l + P 0 

T 2P 0 

(F9) 


The axial diquark case —that we have omitted for 
simplicity— is straightforward. This case differs in the 
explicit appearance of Dirac indices in Eq. (67). In ad¬ 
dition, note that Eq. (67) carries an opposite sign to 
Eq. (69), but this cancels with the different sign in the 
diquark propagators [cf. Eqs. (52 and (54)). We have ne¬ 
glected the term proportional to p^p u in Eq. (54), which 
is suppressed by the diquark mass squared. 


and the Fermi and Bose functions np(E k ) 

( e ^/T + l)- 1 , nB (e k ) = { e ^/T-iy 1 . 
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